ANÁLISIS ESTADÍSTICO PRODUCCIÓN MAÍZ

3.1 Carga de Librerías y Segmentación del Cultivo

Nuestra historia comienza con un universo de información sobre toda la producción agrícola de la empresa. Para encontrar las respuestas sobre el maíz, nuestro primer paso es enfocar el lente, aislando la señal de nuestro cultivo de interés del ruido generado por los demás.

Explicación del código:Se cargan las librerías que nos acompañarán en este viaje analítico. Luego, se leen los datos y, mediante filter(), se segmenta el dataframe para crear una nueva base de datos que contiene exclusivamente los registros de MAÍZ

library(dplyr);library(stringr);library(ggplot2);library(tidyverse);
library(patchwork);library(tidyr);library(naniar);library(corrplot)
library(mice);library(plotly);library(readxl);library(table1)

datos <- readxl::read_excel("C:/Users/johnn/Desktop/GESTION DE DATOS/TALLER PREPROCESAMIENTO/DATOS/eva_df_2025.xlsx")
level_CULTIVO <- c(MAIZ = "MAIZ", maiz = "MAIZ")
datos <- datos %>% mutate(CULTIVO = recode(CULTIVO, !!!level_CULTIVO))
datos <- datos %>% filter(CULTIVO == "MAIZ")

3.2 Carga de Librerías y Segmentación del Cultivo

Los nombres de las columnas originales son poco intuitivos, lo que puede llevar a errores. Para trabajar de manera eficiente y clara, debemos establecer un lenguaje común.

Explicación del código: Se define un vector con nombres claros y estandarizados, y se asigna directamente a las columnas del dataframe datos. Ahora, todos en el equipo sabemos exactamente a qué se refiere cada variable.

nuevos_nombres <- c("departamento", "municipio", "grupo", "cultivo", "año",
                    "area_sembrada", "area_cosechada", "t_produccion",
                    "estado_fisico", "ciclo_cultivo")
names(datos) <- nuevos_nombres

3.3 Creación de la Variable Rendimiento

Producir mucho no es sinónimo de ser eficiente. La empresa necesita una métrica que mida la productividad por hectárea. Aquí es donde creamos al protagonista de nuestro análisis: la variable rendimiento.

Explicación del código: Usando ‘mutate()’, creamos la columna rendimiento. La fórmula t_produccion / area_cosechada se implementa con case_when() para manejar de forma robusta casos especiales, como divisiones por cero o datos faltantes, asegurando la calidad de nuestro indicador clave.

  datos <- datos %>%
  mutate(
    rendimiento = case_when(
      is.na(t_produccion) | is.na(area_cosechada) ~ NA_real_,
      area_cosechada == 0 ~ 0,
      TRUE ~ t_produccion / area_cosechada
    )
  )

3.4 Imputación de Datos Faltantes

Para asegurar la calidad y robustez del análisis, se realizó un preprocesamiento metódico de los datos. Primero, se estandarizaron los nombres de los departamentos a mayúsculas para unificar los registros. A continuación, para tratar los datos faltantes, se evaluaron tres estrategias: ignorar los registros, imputar con la mediana y utilizar una imputación múltiple con el método MICE.

# 1. Estandarizar texto a mayúsculas
# Asumimos que partimos de 'datos_solo_maiz' y lo llamamos 'datos'
datos <- datos %>% 
  mutate(departamento = str_to_upper(departamento))

# 2. Imputación de datos faltantes con MICE
imputR <- mice(datos, m = 5, maxit = 5, seed = 123, print = FALSE) 
Datos_ImputR <- complete(imputR)

# Verificamos que ya no hay datos faltantes
sum(is.na(Datos_ImputR))

## [1] 0

Análisis de los Resultados Se observó que la imputación por la mediana, si bien completaba los datos, tendía a distorsionar la distribución al elevar el impacto de los outliers. Por otro lado, la imputación con MICE generó valores que no producían estimaciones negativas y cuyo comportamiento global era muy similar al de simplemente ignorar los datos. Dado que MICE nos permite conservar la totalidad de los registros sin introducir sesgos significativos, se eligió como el método definitivo. El resultado es Datos_ImputR, nuestro set de datos final: completo, limpio y confiable para el análisis.

4.1 Distribución de la muestra por departamentos

Antes de analizar el rendimiento, necesitamos mapear el terreno: ¿de dónde proviene nuestra información? Este mapa nos muestra la concentración de datos a nivel nacional.

#obtencion de frecuencia en la muestra
datos_freq <- Datos_ImputR %>%
  count(departamento)

#grafico de distribucion de muestra
plot_ly(datos_freq, labels = ~departamento, values = ~n, type = "pie") %>%
  layout(title = list(text = "<b>DISTRIBUCIÓN DE LA MUESTRA POR DEPARTAMENTOS</b>",
      font = list(color = "darkred",  size = 15),x = 0.5, xanchor = "center"))

Análisis de los Resultados La distribución de la muestra está bien repartida entre los departamentos, lo que es positivo para el análisis, ya que las conclusiones generales no estarán fuertemente sesgadas por la información de una sola región.

Llama la atención que Valle del Cauca, a pesar de ser el departamento con la menor cantidad de registros en este conjunto de datos, fue el que demostró tener el rendimiento más alto y consistente en los análisis anteriores. Esto sugiere que su alta productividad no es un artefacto de tener una muestra más grande, sino una característica real de la región.

4.2 Comparación del rendimiento en 2007 vs. 2017

El negocio no es estático. Para entender el presente, debemos mirar al pasado. Este análisis nos cuenta la historia de una década de cambio.

Explicación del código: Filtramos los datos para los años 2007 y 2017. geom_col con position = “dodge” crea un gráfico de barras comparativo que muestra la evolución del rendimiento mediano en cada departamento.

datos_comparacion <- Datos_ImputR %>%
  filter(año %in% c(2007, 2017)) %>%
  group_by(departamento, año) %>%
  summarise(rendimiento_mediano = median(rendimiento, na.rm = TRUE), .groups = 'drop')

ggplot(datos_comparacion, aes(x = rendimiento_mediano, y = reorder(departamento, rendimiento_mediano), fill = as.factor(año))) +
  geom_col(position = "dodge") +
  scale_fill_manual(name = "Año", values = c("2007" = "orange", "2017" = "skyblue")) +
  labs(title = "3. Comparación de Rendimiento Mediano: 2007 vs. 2017", x = "Rendimiento Mediano (Ton/Ha)", y = "Departamento") +
  theme_minimal()

Análisis de los Resultados [Espacio para su análisis e interpretación. ¿Qué departamentos mejoraron su rendimiento en esa década? ¿Alguno empeoró? ¿Qué cambios importantes se observan?]

4.3 Tendencia del rendimiento a través de los años

Esta visualización nos cuenta la historia completa a través del tiempo. ¿Qué departamentos han mostrado un crecimiento sostenido? ¿Cuáles se han estancado o han sido erráticos?

Explicación del código: Agrupamos por año y departamento para calcular el rendimiento mediano anual. geom_line conecta los puntos a lo largo del tiempo, revelando la trayectoria de cada departamento.

tendencia_anual <- Datos_ImputR %>%
  group_by(año, departamento) %>%
  summarise(rendimiento_mediano = median(rendimiento, na.rm = TRUE), .groups = 'drop')

ggplot(tendencia_anual, aes(x = año, y = rendimiento_mediano, color = departamento)) +
  geom_line(linewidth = 1) +
  geom_point(size = 2) +
  scale_x_continuous(breaks = seq(min(tendencia_anual$año), max(tendencia_anual$año), by = 2)) +
  labs(title = "4. Tendencia Anual del Rendimiento por Departamento", x = "Año", y = "Rendimiento Mediano (Ton/Ha)") +
  theme_minimal() +
  theme(legend.position = "bottom")

Análisis de los Resultados [Espacio para su análisis e interpretación. ¿Qué departamentos mejoraron su rendimiento en esa década? ¿Alguno empeoró? ¿Qué cambios importantes se observan?]

4.4 Tablero Gráfico Resumen

Para la toma de decisiones estratégicas, los directivos necesitan una visión global y rápida. Este tablero de mando consolida nuestros hallazgos más importantes en una sola vista.

Explicación del código: La librería patchwork nos permite componer un tablero combinando nuestros gráficos (g1, g2, g4) con una sintaxis simple e intuitiva (+ y /), como si armáramos un rompecabezas.

g1_dash <- Datos_ImputR %>% count(departamento, sort = TRUE) %>% mutate(departamento = reorder(departamento, n)) %>% ggplot(aes(x = n, y = departamento)) + geom_col(fill = "steelblue") + labs(subtitle = "1. Registros por Departamento", x = NULL, y = NULL) + theme_minimal()
g2_dash <- ggplot(Datos_ImputR, aes(x = reorder(departamento, rendimiento, FUN = median), y = rendimiento)) + geom_boxplot(fill = "skyblue") + coord_flip() + labs(subtitle = "2. Distribución del Rendimiento", x = NULL, y = NULL) + theme_minimal()
g4_dash <- ggplot(tendencia_anual, aes(x = año, y = rendimiento_mediano, color = departamento)) + geom_line(linewidth = 1) + labs(subtitle = "4. Tendencia Anual del Rendimiento", x = "Año", y = "Rendimiento (Ton/Ha)") + theme_minimal() + theme(legend.position = "none")
(g1_dash + g2_dash) / g4_dash + plot_annotation(title = 'Tablero Gráfico Resumen del Cultivo de Maíz')

4.5 Tabla de Indicadores Descriptivos

Mientras que los gráficos nos dan la intuición, las decisiones de negocio requieren los números detrás de las imágenes. Esta tabla proporciona las estadísticas descriptivas precisas.

Explicación del código: table1() es una potente función que genera tablas con calidad de publicación. Con una simple fórmula (~ rendimiento | departamento), obtenemos un resumen estadístico completo, desglosado por departamento.

table1(~ rendimiento | departamento, data = Datos_ImputR, topclass="Rtable1-zebra")

Análisis de los Resultados [Espacio para su análisis e interpretación de la tabla. Puede comparar directamente la media y la desviación estándar entre departamentos, o comentar sobre los rangos (mínimo y máximo) para reforzar las conclusiones de los boxplots.]